LINEAR PROGRAMMING BOUNDS FOR CODES IN 
GRASSMANNIAN SPACES 



CHRISTINE BACHOC 

Abstract. We develop the linear programming method in order to 
obtain bounds for the cardinality of Grassmannian codes endowed with 
the chordal distance. We obtain a bound and its asymptotic version that 
generalize the well-known bound for codes in the real projective space 
obtained by Kabatyanskiy and Levenshtein, and improve the Hamming 
bound for sufficiently large minimal distances. 



1. Introduction 

Philippe Delsarte has introduced the so-called linear programming method, 
in order to find bounds for the size of codes with prescribed minimal dis- 
tance, in the classical case of codes over finite fields. This method, also 
called Delsarte method or polynomial method, exploits a certain family of 
orthogonal polynomials attached to the situation, the Krawtchouk polyno- 
mials, and their positivity property. These polynomials and their properties 
are intimately related to the action of the symmetric group on the Hamming 
space. Delsarte method has proved to be very powerful, and was extended to 
many other situations, where the underlying space is symmetric of rank one, 
and is homogeneous under the action of a certain group of transformations. 
Examples of such spaces are: the Johnson space, the Grassmannian space 
over a finite field, the unit sphere of the Euclidean space, the projective 
spaces over the real, complex and quaternionic fields (for these last spaces 
see P, 0). 

In recent years, codes over the real Grassmannian space have attracted 
attention, motivated by their application to information theory, more pre- 
cisely to the so-called space-time codes, used for multi-antenna systems of 
communication. The distance usually considered is the chordal distance, 
introduced in [3.. , and defined in the following way (more details are given in 
the next subsection): The Grassmannian space of m-dimensional subspaces 
of R n , where m < n/2, is denoted by Q m ,n] to a pair (p,q) of elements of 
G m ,n is associated m principal angles 6\, . . . , 8 m £ [0, vr/2]. Let in := cos 2 0j. 
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d c (p,q) ■= , ^ sin 2 Oj 



In 0, the authors give bounds for the size of Grassmannian codes, called 
the simplex and orthoplex bounds. The main drawback of these bounds is 
that they are only valid in a certain range of minimal distances. In [I], an 
asymptotic bound, derived from the Hamming bound, is given. Another ap- 
proach is developed in 0, where bounds are given for codes whose principal 
angles are subject to certain constraints (the so-called /-codes), which arise 
naturally from the notion of Grassmannian designs introduced in [H]. 

In this paper, we extend Delsarte method to the Grassmannian codes, 
exploiting the zonal polynomials attached to Q m ,n- These are symmet- 
ric polynomials in the m variables yi, . . . ,y m ; they belong to the family 
of orthogonal generalized Jacobi polynomials (see the next subsection). In 
the second section, we recall, or settle the properties of these polynomials 
needed to perform linear programming bounds; these properties are easy 
to obtain by straightforward generalization of the arguments used in the 
classical cases. In fact, the principles underlying the LP method would re- 
main true for the zonal polynomials attached to any symmetric space. The 
real difficulties start when one wants to actually perform explicit bounds, 
because the polynomials have (for m > 2) several variables. The low degree 
cases are still easy to manage; this is done in section |H1 where we recover the 
simplex bound as the bound arising from the case of degree one, and give 
new bounds from polynomials of degree 2 and 3. In the forth section, we 
propose a strategy based on the eigenvalues of certain symmetric endomor- 
phisms, which extends the one variable method based on the zeros of the 
polynomials and on Christoffel-Darboux formula, but avoids to deal with 
zeros of polynomials in several variables. We obtain an upper bound for the 
size of a code C with minimal distance 5, which is expressed in terms of 
the largest eigenvalue fTheorem 14.41 and Corollary 0) . Section \5\ settles the 
asymptotic behavior of this largest eigenvalue (Theorem l5.3|) . and in section 
IHlwe derive the following asymptotic version of the bound: 

Theorem 1.1. Let C be a code in Qm n with minimal chordal distance S , 
let s := m — S 2 G]0, m[ and let 




Then, when n — > +oo 





< m((l + p) log(l +p) - plog(p)). 



Our bound coincides with the bound given by G. Kabatiansky and V. 
Levenshtein in [7] for the case of the real projective space, corresponding 
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to m = 1. But it beats the Hamming bound of [1] only when the minimal 
distance is relatively big. 

1.1. Basic facts about Grassmannian spaces and their zonal poly- 
nomials. We repeat here, without proofs, some well-known facts about 
Grassmannian spaces and their zonal polynomials. Some useful references 
for the mathematical background are: |Hj, [2] for the representations of the 
orthogonal group, ^3], |15j . for the Grassmannian spaces and harmonic 
analysis on it, [HI- [HI f° r multivariate orthogonal polynomials. 

The real Grassmannian space, denoted by Q mn (m < re/2), is the set 
of rre-dimensional R- linear subspaces of R n . The orthogonal group 0(re, R) 
acts transitively on Q m ,n] a transformation stabilizing a given element po 
also stabilizes its orthogonal complement p$ and therefore the stabilizer of 
Po is isomorphic to the direct product 0(m,R) x 0(re — m,R). Hence we 
derive the identification of Q m ,n with the set of classes: 

Q m)n ~ 0(n,R)/(0(m,R) x 0(n-m,R)) 

from which Q m ,n inherits the structure of a (compact) differential variety, 
and a 0(re, R)-invariant measure that will be normalized so that Jg dp = 
1. It is worth noticing that the case m = 1 corresponds to the real projective 
space. 

In order to understand the action of 0(re,R) on pairs (p,q) S G m ,n 2 , we 
need to introduce the principal angles between p and q. These are m angles 
9\, . . . , 6 m £ [0, 7r/2] defined in the following way: 

Let p\ C p, qi C q be two lines such that the angle 0\ between p\ and q\ 
is minimal. If m = 1 we have finished, otherwise let p' be the orthogonal 
complement of p\ in p, q' be the orthogonal complement of q\ in q; we define 
recursively 62, ■ ■ ■ , 9 m to be the principal angles associated to the pair (p' , q') 
in G m -i,n- We introduce the notation := cos 2 ^; when needed, we may 
denote rather yi(p,q), 9i(p,q). A classical result on the geometry of Q m ,n is 
the following: 

Proposition 1.2. Two pairs (p,q) and (p',q') are in the same orbit under 
the action of the group 0(re,R), i.e. there exists a £ 0(re,R) such that 
o~(p) = p' and o~{q) = q' , if and only if 

yi(p, q) = yi(p', q') for all 1 < i < m. 

The previous proposition expresses the fact that the orbits under the 
action of 0(re, R) of the pairs (p, q) G G m ,n 2 are characterized by the m- 
tuple of real numbers (yi(p, q), ■ ■ ■ , y m (p, q))- It becomes clear that, for 
rei > 2, Q m ,n is not 2-point homogeneous, i.e. a single distance on Q m ,n 
cannot characterize these orbits (while it is the case for other spaces of 
interest in coding theory, like the Hamming and binary Johnson spaces, or 
the unit sphere of the Euclidean space). It is the reason why we shall deal 
with mutivariate polynomials. Also, it shows that the choice of a distance 
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on G m ,n is sort of arbitrary. We shall stick to the chordal distance in this 
paper, as introduced in [H]: 



d c (p,q) ■-- 



sin 2 ■ 



n£ ■ \ 



E: 



m 

i=l 



Other possibilities are the Riemannian distance v/2^™i#f which be- 
haves somewhat badly because it is not smooth; the max distance maxj#j, 
etc.. The "product distance" (which is not a distance in the metric sense) 
( Yii sin 8j) seems to be relevant in the context of space time codes. 

Now we consider the space L 2 {Q m ^ n ) of functions / : G m , n —> C such that 
Jg \f\ 2 dp < +00. This is a C-vector space, endowed with the hermitian 
product: 



<f,g>= / f{p)g{p)dp 



and with the left action of the orthogonal group given by: 

(vf)(p) = f(°-\p)) 

(for which the above hermitian product is of course invariant). 

Its associated algebra of zonal functions (also called the Hecke algebra) 

is: 

Z := {Z :Q m ,n -> C | Z(p, •), Z(-,q) £ L 2 {Q m ^ n ) and 
Z(a(p),a(q)) = Z(p,q) for all a G 0(n,R)} 

Form Proposition 11.21 since Z £ Z is constant on the orbits of 0(n,M) on 
G m ,n 2 , it can be given the form: Z(p, q) = z(y 1 (p, q), . . . , y m (p, q)) for some 
function z. 

The explicit decomposition into 0(n, M)-irreducible subspaces of 
L 2 (Gm,n), and the corresponding structure of Z, where investigated for the 
first time by James and Constantine (^1]). It is now a standard result on 
the representation of the classical groups (see [0]). 

Recall that the irreducible representations of 0(n, M.) are (up to a power 
of the determinant) naturally indexed by partitions k = . . . , n n ), where 
k,i > ■ ■ ■ > n n > (we may omit the last parts if they are equal to 0) . 
Following [Hj, let them be denoted by V£. For example, Vn = CI, and 
= Harmfc the space of homogeneous of degree k, harmonic polynomials 
in n variables. 

The length £(k) of a partition k is the number of its non zero parts, and 
its degree deg(fc) also denoted by equals Yl7=l Ki - 
Then, the decomposition of L 2 {Q min ) is as follows: 



L\G m ,n) - ©K 



2k 
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where k runs over the partitions of length at most m and Ik stands for 
(2rei, . . . , 2K m ), meaning that only partitions with even parts enter the de- 
composition. We can see that the multiplicities in this decomposition are 
all equal to one, which translates the fact that the space G m ,n is a sym- 
metric space. Consequently, to each irreducible component V^ K is asso- 
ciated a uniquely determined (up to a normalizing factor) zonal function 
P K (yi, . . . , y m ), in the sense that 

p (resp. q) i-> P K (yi{p, <?),..., y m (p, <?)) G V^ K 

and 

It turns out that the P K are symmetric polynomials in the m variables 
... j y m , of degree |«| 3 with rational coefficients once they are normalized 
by the condition P K (l, . . . , 1) = 1. Moreover, the set (P K )\ K \<k 1S a basis of 
the space of symmetric polynomials in the variables y±, ... , y m of degree at 
most equal to k, denoted by 

Since the irreducible subspaces of L 2 {Q m ^ n ) are pairwise non isomorphic, 
they are orthogonal for the 0(n, M)-invariant hermitian product defined 
above. This hermitian product induces an hermitian product on the space 
of symmetric polynomials, denoted by [, ] , for which the polynomials P K are 
orthogonal. More precisely, it is given by the positive measure, calculated 
in [TU, 

m m 

dfi = A J] \Vi - Vi\ JI%" 1/2 (! - yiT /2 ~ m ~ 1/2 dyi 

i,j=l 1=1 
i<j 

(where A is chosen so that Jj ^ m dfi(y) = 1). and 

[/,<?]= / f(y)~gTy)dn(y). 

J[0,l] m 

One recognizes a special case of the orthogonal measure associated to 
generalized Jacobi polynomials 

We let Ilfc be the subspace of generated by the polynomials P K , with 
\k\ = k, so that we have the orthogonal decomposition: 

Sk = Sk-i J- Ilfc. 

Let the dimensions of Sk, n& be denoted respectively by Sk, vr^. The number 
7Tfc is also equal to the number of partitions k of k in at most m parts. These 
dimensions also depend on m, although it does not reflect on our notation, 
for the sake of simplicity. 

In view of the explicit calculation of the polynomials P K , it is better to 
use the following characterization, which involves the polynomials C K , which 
are themselves the zonal polynomials associated to the symmetric space 
GL(m, R)/0(m,R) (these polynomials are Jack polynomials, normalized by 
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C K (1,...,1) = 1, see [H], ^21)) and the differential operator A induced on 
C[yi, . . . ,y m ] by the Laplace Beltrami operator of Gm,n- The condition: 
for all 1 < % < m, m > \x% is denoted by: k> fi. 

(i) P K is an eigenvector for the operator 



m Q2 m Q 

4 = 1 1 l^] = l 

m ^ m rv) 

+ (- -m + 1) > yi- >^yi^ 

2 ^ ^ ! ^ 

m d 1 m d 

- E «.(»- ar'ajj-jE as 

l?frj=l 1=1 

(ii) P K = /3 K C K + X^i/o/i Pk,ijPh 

(iii) P«(l, . . . , 1) = 1. 

Condition (ii) is needed to avoid the multiplicities of the operator A. 

Examples: the effective computation of the polynomials P K following the 
method described above leads to, up to the normalization imposed by (iii): 



P = l 

\l) = si 



„ m 
Pi . 

n 

\2 ™2(^_i\2 



(m — 1) m (m — 1) 

P (H) = ai ~ 81 + 



n-2 2(n-l)(n-2) 
2 2(m + 2) 2 m 2 (m + 2) 2 



^(2) =S2 + -<?1 ~ Sl + 



3(n + 4) 3(n + 2)(n + 4) 
where si = £i<i< ro W> s 2 = Ei<i< m % ? > 

°"i = Yli<i<j< m yiyj- 

Remark: The complex Grassmannian £/ TOin (C) is more commonly used 
in the context of space-time coding. It affords the transitive action of the 
unitary group U(C, n); similarly one defines principal angles (0\, . . . , 9 m ) be- 
tween two elements of Gm,n(C). The U(C,n) decomposition of L 2 {Q m ^ n {<C)) 
and the associated zonal polynomials are computed in Jl] so one can play 
the same game concerning bounds of codes. On the other hand, a bound 
is also obtained from the embedding G m< n(C) C G2m,2n(^) (if (#i> • • • ,0 m ) 
are the principal angles associated to a pair (p,q) of elements in Gm,n(^), 
the 2m principal angles associated to the pair (p, q) , seen as elements of 
^2m,2n (M) , are simply (0 X , 0\ , 2 , 2 ,-••)) • 
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2. Zonal polynomials associated to Q m ,n and the LP bound 

In this section, we settle the properties of the polynomials P K relevant 
for the LP bound, settle this bound, and show how the Christoffel-Darboux 
formula can be exploited in that context. 

The dimension of V£ is denoted by d K . Explicit formulas for d K can be 
found in |S]; however we do not need them before Section 5. 

Proposition 2.1. The polynomials P K , normalized by the condition 
P K (1, . . . , 1) = I, satisfy: 

(i) [P K ,P K ]=d£ 

(ii) (Positivity property): For all finite set C C Qm,n, 

p,q&C 

(hi) Let p v K lL be defined by the property: 

PkP^i = Pn^Py 

V 

The numbers p v K „ are non-negative numbers. 

Proof. These properties where already pointed out in [5J [Lemma 2.2] and 
step on very general arguments (see {TBI [Theorem 3.1]). For the sake of 
completeness, we briefly recall the arguments. Let e%, . . . ,ed 2K be any or- 
thonormal basis of the subspace H^ n of L 2 (G m ,n) isomorphic to V^ K . Let 
P K (p, q) := P K (yi(p, q), . . . , y m (p, <?))• It is well known that we have (this is 
called the addition formula) 

^ din 

Pk(p,q) = -i— y^ei(p)ei(g). 

i=l 

As a consequence, from the expression 

[P K ,P*] = [ P K (p,q)P,(q,p)dq 
(i) follows. Moreover, 



p,q£C p,q€C i=l 

-, d2 K 
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hence (ii). More generally, for any function a : C — > C, we have: 
oc(p)a(q)P K (p, q) =— I a (P) e i(P) 

p,</GC 2K i=l p£C 

> 0. 

Conversely, assume F G S k is a polynomial with real coefficients, such that, 
for any finite set C C G m ,n and any function a : C — > C, 

£ a{p)c4q)F(p,q) > 0, 

and let us prove that F expands on the P K with non-negative coefficients. 
Taking limits, we have, for any a G L 2 (Q m ^ n ), 



a(p)a(q)F(p, q)dpdq > 

and hence, using the addition formula, 

P K (p,q)F(p,q)dpdq > 0. 



Qm,n 

If F = J2\ u \<k fvPvi the left hand-side equals f K /d2 K , which proves that 
the coefficients f K are non-negative numbers. Using once again the addition 
formula, it is easy to show that the product P^P^ holds this general positivity 
property, and therefore expands on the P v with non-negative coefficients. 

2.1. The principles of the LP bound. The positivity property of the 
polynomials P K is the basis of the linear programming method to upper 
bound the cardinality of <5-codes. 

Definition 2.2. A Grassmannian code C satisfying the constraint: 
For allp^ q G C 2 , d c (p, q) > 5. 

is called a 5-code. 



Proposition 2.3. Assume F k £ satisfy: 

(i) F k = J2\ K \<k f* P K with fn>0 for all k, f > 

(ii) F k (yi, ■ ■ ■ ,y m ) < for all (y u . . . ,y m ) € [0, l] m such that 

YT=i Vi<m-5 2 

Then, the following bound holds for the cardinality \C\ of any 5-code: 

JO 
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Proof. This is a standard argument, that we recall here. Let C be a <5-code. 
As before, we let F k (p, q) = F k (yx(p, q), . . . , y m (p, q)). We calculate 

£ F k (p,q)=J2 A( E 

Assumption (ii) leads to F k (p,q) < when p ^ q. The remaining terms 
of the left hand-side, corresponding to p = q, contribute by |C|i^(l, . . . , 1). 
Assumption (i), together with the positivity property of the polynomials 
P K (Proposition 12. ll (ii)), show that all the terms of the right hand-side are 
non-negative. When k = (0), P K = 1 and the contribution is /o|C| 2 . We 
obtain 

|C7|F fc (l,...,l)>/ |C7| 2 

equivalently 

JO 

It is worth noticing that equality in this inequality happens if and only if, for 
all 1 < \k\ < k such that f K ^ 0, V, qeC P K (p, q) = and, for all p ^ q € C, 

-^fc(P) q) = 0- The first condition says that C is a 2fc-design in the sense of [£] 
(when it holds for all 1 < \k\ < k), and the second one that C is an i^-code 
in the sense of 0. 

2.2. The three-term relation and the Christoffel-Darboux formula. 

We join here more material on the sequence of polynomials P K , that will be 
of later use. The results presented here are essentially established in ex- 
cept that we deal with symmetric polynomials. Following [13] . the (column) 
vector of the polynomials P K with \k\ = k is denoted by ¥ k . If necessary, we 
order the partitions of the same degree in increasing lexicographic order. 
We also set 

o- ■= yi + 2/2 H Vy m 

and, when necessary, we make the involved variables explicit, by writing 
o~(y) rather than a. The ir^ x ix k diagonal matrix, denoted by D k , with 
entries 

D k [ K ,K] = d 2K := dim(V* K ) 

is the inverse of the Gram matrix of ¥ k . 

Next result is an analogue of the so-called "three-term relation" . 

Theorem 2.4. For all k > 1, there exists matrices A k , B k , C k , of size 
respectively ir k x TT k +i, ir k x Tr k , ir k x iik-i, such that: 

aF k = A k F k+1 + B k F k + C fc P fc _i. 

Moreover, {D k B k y = D k B k and D k C k = (Z) fe _iA fc _i) t . 



10 



CHRISTINE BACHOC 



Proof. The polynomials aP K with \k\ = k are symmetric of total degree k+1 
so they afford a decomposition over the (^)|d<fe+i- Moreover, [aP K ,P^\ = 
[P K ,aP M ] = if \fjt\ < k- 2. 

If = |/c| = k, we have: £*.[«, //] [P M , P M ] = [o-P K ,P M ] = [P K ,oP M ] = 
Bi.[fj,, k][P k , P k ], which proves that the matrix D^B^ is symmetric. The 
same argument shows that D^C^ = (P^-i^-i)* • 

Notations: We want to define «W (respectively /%)) to be the partition 
obtained from k by increasing (respectively decreasing) the i-th. part «j by 
one. This is not possible for all i, since the result should be also a partition, 
i.e. the new parts should be in decreasing order. Hence we define (where 
K m +i := 0): 

u(k) := {1} U {i £ [2, m] \ Ki-i > Ki} 
d{n) := {i G [l,m] | m > 

The set u(k) is the set of indices i for which «W makes sense (respectively 
for Moreover, if |/t| = k, \k®\ = k + 1 and \ku\ \ = k — 1. 

Otherwise explicitly mentioned, in the rest of this paper, k, k! are par- 
titions of degree k, while /x, // are partitions of degree k + 1 and u, u are 
partitions of degree k — 1. 

Proposition 2.5. T/ie following properties hold: 

(i) Por a// k, fx and k' , Aj~ [n, fj] > ; and P^f/t, k'] > 0. 

(ii) T/ie coefficients of the matrix are equal to zero, except the coef- 
ficients Ak[ W], which are positive. 



Proof. The first assertions are equivalent to: [oP^PJ > for all k, /U of 
any degree. But [aP^P^] = [l,aP K P M ] and a = m(l - S)Pi + Joint 
with Proposition 12. If iii) ■ we obtain [l,oP K P^] > 0. 

The coefficients Ak[n,(j] can be more precisely calculated, using ^ [T71 
Lemma 7.5.7]). Since we do not normalize the polynomials C K in the same 
way, we introduce coemcients M such that 



They differ by a positive multiplicative factor from the generalized binomial 
coemcients (^) defined in jTJj; see also Then we have 



A k [K,Ll] 



It is known that the generalized binomial coefficients 0£) are equal to zero 
when // is not equal to one of the «W; consequently the same holds for 



A fc [/c,/i]. Moreover, since (^°) > 0, also [ K ^'} > and A fe [re,«W] / 0. 
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Theorem 2.6 (Christoffel-Darboux Formula). Let 

Qk-= ^2 A k[K,l4 P v e n fc+l- 

\^\=k+l 

With the previous notations, we have: 

(i) For all k > 0, 

£ d 2v P v {x)P v {y) = 
\v\<k 

a(x) - a{y) 

(ii) Moreover, if e := YT=i ^> 

\u\<k 

£ ^([eQ K (y))P K (y) - [eP K (y))Q K (y) 

\n\=k 



Proof. The proof of (i) is the same as |13| [Theorem 3.5.3]. Note that we 
cannot hope for a formula for each X{ like in |13j . since we should stick to 
symmetric polynomials. If 

S s : = (A s P s+1 (aO)'£yP s (y) - P s (x)' AAP s +i(y) 
= d 2n{Q,{x)P K {y) -P K (x)Q K (y)), 

\k\=s 

from the "three-term relation" of Theorem 12.41 we have: 

S s - S s _! = (o-(x) - <x(y))P s (aO'A,P a (y). 
The formula (i) follows from summing up these identities, for 1 < s < k. 
In the equation (i), we replace Q K {x)P K {y) - Pn(x)Q K (y) by 

Q K (x)(P K (y) - P K (x)) - P K (x)(Q K (y) - Q K (x)). 
Then, if we specialize x 2 = y 2 , ■ ■ ■ , %m = Vm and let x\ tend to yi, we obtain 

^ d 2 „(P„(y)) 2 = 

\u\<k 

£ d 2K A k [K,fA({^f(y))P K (y) - {^{y))QM 

\n\=k 

|/Lt|=fe+1 
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The same identity holds when one replaces y\ by any yj ; if we sum up all 
these identities, we obtain the more symmetric formula (ii). 

Remark 2.7. The left hand side of the Christoffel-Darboux formula 

Kk(x,y) := ^2 d 2v P v (x)P v (y). 
\u\<k 

is the reproducing kernel of the space of symmetric polynomials of degree at 
most k. It satisfies the characteristic property: for all Q £ S k , [K k (x, .), Q] = 
Q(x). 

2.3. An LP bound from Christoffel-Darboux formula. In the classical 
cases, Christoffel-Darboux formula is involved in the setting up of bounds 
of the type \C\ < M(5) where M{5) is an explicit function of 5. Usually 
the running interval of 5 is divided into subintervals, related to the zeros of 
the zonal polynomials. This is the line followed in ^H], and also in [7]; see 
|19j for a unified presentation. In this section, we follow this method, and 
analyze the difficulties arising from the several variables situation. 

The numerator, of degree k + 1, of the right hand side of Christoffel- 
Darboux formula (Theorem I2.6f i^ is denoted by N k+ i(x,y). We consider 
the polynomial in the variables yi, . . . ,y m , of degree 2k + 1, 

F2k+i{x,y) : = -N k+1 (x,y)K k (x,y) 
= N k+1 (x,y) 2 
a(y) -a{x)' 

In order to make clear that only the yx,...,y m are variables, while the 
will specialize to real values, we denote it by F 2k +i(x, •)• 

Proposition 2.8. Let x £ [0, l] m satisfy a(x) > s := m — 5 2 . Assume the 
following conditions hold: 

(i) For all k, \k\ < k, P K (x) > 

(ii) For all k, \k\ = k, Q K (x) < 

Then, F2 k +i{x,-) satisfies the conditions required in Proposition^^ 



Proof. We have: 

N k+ i{x,yf 



F2k+i(x,y) 



a(y) - a(x) 

hence condition (ii) is satisfied when s < o~(x). 

To prove condition (i), we point out that, if F and G are two polynomials 
with non- negative coefficients on the P K , then the product FG holds the 
same property. This is a direct consequence of Proposition 12. in n). 
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From the definition of K k (x,y), its coefficient on P v with \v\ < k equals 
d 2u P v (x) (and for higher degree partitions it is zero). On the other hand, 

-N k+1 (x,y) = d 2K {P K (x)Q K (y) - Q K (x)P K (y)) 

\n\=k 

|/i|=fc+l |k|=/c 

- d 2K Q K (x)P K (y) 

\n\=k 

The coefficient A k [K, fj] is always non-negative. Clearly, under the condi- 
tions of the proposition, the coefficients of — N k+ i(x,y) on the P K and P^ 
are non-negative. 



Corollary 2.9. Assume x satisfies the conditions of Proposition ^^ . Then, 
for all S-code C, 

I ll) — (J{ x i ] ( > i. , <!■■>,, /'„(.)•)'" 

\c\< 



(m - a(x)) ( J2\u\<k d2uPu(x)Y 



J2\ K \=k d2KPn{x)Q K {x) 



Proof. In order to apply Proposition 12.31 we are left with the computation 
of f and of F 2k +i(x, (1, . . . , 1)). Since F 2k +i(x,y) = K k (x,y) 2 (a(y)-a(x)), 
we have 

F 2k+1 (x, (1, . . . , 1)) = ( d 2u P u (x))\m - a(x)). 

\u\<k 

Using the orthogonality of the P K , we obtain 

fo ■ = [F 2k +i(x, •), 1] = ~[K k (x, -),N k+1 (x, •)] 

= -[J2 d 2u Pu{x)Pv, d 2,Qn(x)P K ] 
\v\<k \n\=k 

= ~ d 2K P K (x)Q K (x). 

|fc|=fc 

The main problem with this approach, is that, in general, we don't even 
know if the inequalities (i) and (ii) of Proposition 12.81 have a solution x. 
In case these inequalities define a non empty area of HL m , a second problem 
would be to optimize the choice of x in this area. In the classical case m = 1 , 
Qk = Pk+i (up to a positive multiplicative factor). The interlacing property 
of the real zeros of the orthogonal polynomials P k , ensures that one can take 
x G [z kl z k+ i], where z k is the largest zero of P k , so that P k+ i(x) < and 
Pi(x) > for all i < k. Moreover, one uses asymptotic estimates of these 
zeros to derive an asymptotic bound for the size of codes. 

In the general case m > 2, we don't have such tools to deal with the 
inequalities of Proposition 12. 81 which seem to be intractable in general. The 
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first case k = 1, leading to a polynomial of degree 3, is however discussed in 
the next section. On the other hand, one can think of the zeros of orthogonal 
polynomials in one variable as being the eigenvalues of the so-called Jacobi 
matrices associated to the sequence of polynomials. We study in section 
0]the eigenvalues of the analogous matrices in the general case, and derive 
bounds for codes, which contain as a special case the bound obtained from 
a possible solution of these inequalities. 



3. LP BOUNDS OF SMALL DEGREE 

We take the following notations: let s := m — 5 2 , the maximal value of 
a among pairs of points of a code C. We are looking for a function M(s) 
such that \C\ < M(s). Obviously, M(s) is an increasing function. In this 
section, we discuss the cases of small degree k, trying to optimize the choice 
of Ff. in Proposition 12.31 

3.1. Degree 1. Let iq = 1 + fiP\, with /i > (condition (i)). We have 

P = - ((T — 177,2 ) 
1 m(n—m) \ n ' ' 

When a £ [0, s], 1 + f\P\ should be non-positive (condition (h)). There- 
fore, The zero of I + /1-P1 should be greater than s. It leads to the condition: 

m? m(n — m) 
s < 



n nf 



Since fi > 0, we obtain the necessary condition s < — . The smallest 
value for f\ is then 

— m(l — m/n) 



corresponding to a polynomial proportional to a — s. We obtain the bound 



m 2 , „, m — s 



if s < — , |C| < 



which is the so-called simplex bound proved in [3]. 

3.2. Degree 2. We restrict ourselves to polynomials which are divisible by 
a — s. Then, such polynomials are polynomials in a. We write: 

F 2 = (a- s)(a -b) = f 2 P 2 + /11P11 + fiPi + /o- 
with the condition that b < 0. With t = s — m 2 /n, we find: 
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h 

hi 



m{m + 2)(n — m)(n — m + 2) 
' 3(n + 2)(n + 4) 

2m{m — l)(n — m)(n — m — 1) 
' 3(n-2)(n- 1) 



/, m\ /m 2 4(n — 2m) 2 

/i :=m 1 — + — i ; - - t-b 

\ nJ \n n(n — 2)(n + 4) 

2m 2 (n — m) 2 m 2 
Jo : =— o? Tw — — * + 

The condition /o > 0, when t > 0, is equivalent to 



. , m 2 ( 2(n — mY 

(2) b > — 1 



n \ tn(n — l)(n + 2) 
(and when i <= is always fulfilled), which implies 



. . 2(n — m) 2 

(3) t < 



n(n - l)(n + 2)' 
The condition fx > is equivalent to 

m 2 4{n — 2m) 2 

(4) b < — + 1 - - t. 

n n(n — 2)(n + 4) 

One can check that the right hand side of (j3J) is positive for m > 2, when 
t satisfies Q. 

The bound B = (f 2 + fn + fi + /o)//o equals 

(5) B = 5 2 ^. 

To 

As a function of b, it is decreasing when t G [ w ^i"(~+2) ' n(n-i)(n+2) [' and 
hence the best choice of b is 6 = 0. We obtain the bound: 



Theorem 3.1. If s £ 



n , 2(ra-m) 2 

U ' n ' n(n-l)(n+2) 



(6) \C\<±' 



m v_ s _|_ m£ _|_ 2(n-m) 2 



n ^ l)(n+2) 

This bound, which is an increasing function of s, improves on the simplex 
bound when s > ^ — ^^^",7!"^ . Their common value at s = ^ 



16 



CHRISTINE BACHOC 



2m(n—m) - q jn+1 
n(n-l)(n+2) 

reads: 



is 



(2 )• However, the orthoplex bound proved in [31 (5.6)] 



s < m /n =>• \C\ < 1 



2 

s = m 2 /n \C\ < (n - l)(n + 2) 

and is better than ® in the range - J(^i){J!+2) > ^ we P m § in 

© the value s = m?/n, we find that \C\ < ^(l-m/n)^ wmc h * s better than 
the orthoplex bound when m < n/2. We recall that the orthoplex bound 
is attained for a family of codes with n = 2 l , m = n/2, constructed in |2U1 
Theorem 1]. These codes are also optimal 6-designs (see |21|). 

3.3. Degree 3. We do not study general polynomials of degree 3 but rather 
apply the approach described in subsection 12.31 

The polynomial F3 has degree 3, and is again a polynomial in a. In the 
following, we calculate the best choice for x (and discuss its existence). Let 
u := cr(x) — m 2 /n. We should have: 

(i) u > s — m 2 /n 

(ii) u > (Condition (i)) 

(iii) « 2 ~ ntt~2)^+A) u ~ n™n-i)(n+2) ^ (Condition (ii)) 

The polynomial of degree 2 occurring in (hi) has a positive discriminant, 
and a unique positive root that we shall denote by U2- Let b and c be the 
coefficients of this polynomial, so that it is equal to u 2 — bu — c, and let 
d : = ^|t^t . The bound is then equal to: 

(n — l)(n + 2)(u + d) 2 (m — u — m 2 /n) 
{U) := 2u{u 2 -bu-c) ' 

The calculation of B'(u) shows that it is increasing in the range [ui,^] 
(the numerator has the form: u + d times a degree 3 polynomial with a 

2 2 

unique real root Ux). Hence, for s 6 [u\ + U2 + ^-], the best choice for 
u is u = s — m 2 /n. 
We obtain: 

n-n o r» r ± 7 2m(n—m) > 4(n-2m) 2 2m 2 (n-m) 2 

Iheorem 3.2. Lei a = —> — Vw — rWr, o = ? ?77 , *\ , c = —t, — ^-ttt — tttt- 

n(n— l)(n+2) ' n(n— 2)(n+4) ' n^(n— l)(n+2) 



//*€ 



m 2 m 2 6 /fe 2 
ra n 2 V 4 



|c|< {m-s){s-^+d) 2 (n-l){n + 2) 
| -2(,-^)(-( S -^) 2 + 6( S -^) + c) 
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4. The endomorphisms T k 

We introduce an endomorphism T k : S k — > which eigenvalues will play 
the role of the zeros of the zonal polynomials in the rank one case. 

Proposition 4.1. Let 

Tk ■ S k — > 5^ 

where the orthogonal projection on S k is denoted by pr Sfc (note that, in 
general, aP does not belong to S k but rather to S k +i)- 

The endomorphism Tk is a symmetric endomorphism of Sk, and is an 
isomorphism. 



Proof. We have, for all P,Q G S k , [T k (P),Q] = [aP,Q} = [P,aQ] = 
[P,T k (Q)]. Moreover, [aP,P] = [a,P 2 ] > unless P = 0, because of the 
positivity of the measure on H.[yi, . . . , y m ] Sm ■ Thus Tk is injective. 

Let Jk be the matrix of this endomorphism in the basis {P K , \n\ < k}. 
From the three-term relation (Theorem 12. 41) . Jk is the block-tridiagonal ma- 
trix: 



(7) 



/B A 
Ci B x A x 

C2 B2 A2 



Jk 



\ 



C 3 



V 



/ 



Ck Bk 

It is worth noticing that the matrix J k itself is not symmetric, because the 
polynomials P K are not of norm 1. We shall later introduce and calculate 
the symmetric matrix J' k obtained in the normalized basis. 

In the end, we shall need some very precise information on the coefficients 
of J'k- For the moment, the only, but crucial, property that we will exploit 
is the fact that it is non-negative and irreducible. 



Lemma 4.2. The eigenvalues of Tk are real, and belong to ]0, m[. The 
maximal eigenvalue ofTk, denoted by \ k) is of multiplicity 1, and possesses 
an eigenvector with positive coordinates. Moreover, Ak-i < 



Proof. The matrix Jk is non- negative and irreducible in the sense of |22j . 
because of Proposition ^. 51 (note that the coefficients ^4fc[ft, n^'] are positive). 
Moreover, it is the matrix of a symmetric endomorphism, so its eigenvalues 
are real. From |22[ Perron- Frobenius Theorem], it follows that the maximal 
eigenvalue has multiplicity equal to 1, and that, if v is an eigenvector, either 
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v or — v has positive coordinates. Let us now prove that all its eigenvalues 
belong to ]0, m[. 

For any v G S k , v ^ 0, we have [or>,t>] = J av 2 dfi(y), where d\i is a 
positive measure. We integrate on the domain [0, l] m , on which < a < m, 
hence < [<7v,t>] < m [v,v]. If v is an eigenvector of T k associated with an 
eigenvalue A, we have [av,v] = [Xv,v] = X[v,v], so we can conclude that 
< A < m. 

Now let v be an eigenvector of T k -i for X k -i, assumed to be of norm 1. 
We have 



0~V = Xk-lV + U 

with u G Ilk- Obviously, since deg(av) = 1 + deg(v), v must be of degree 
exactly k — 1 (and u ^ 0). 
Since 

, [T k (x),x] 
Afc = max 



xG5 fc \{0} [X,XJ 

we have [Tk(v),v] < Afc. But [Tfc = [au, «] = Afe_i. The equality 
Afe-i = Afc would mean that f is an eigenvector of Tk, which is not possible 
since it has degree k — 1. 

In the case m = 1, the eigenvalues of are exactly the zeros of the 
polynomial Pfc+i. In the general case, we prove in next lemma that common 
zeros of the polynomials Q K give some of the eigenvalues. However, we do 
not know if such common zeros do exist, neither if all of the eigenvalues are 
obtained that way (and may be it is not so important): 

Lemma 4.3. Let a £ [0, l] m be a common zero of the polynomials 

M=fc+i 

for all k, \k\ = k. Then, v := 5Z|i/|<& d2 V P v (a)P v is an eigenvector ofTj- for 
the eigenvalue a(a). 



Proof. It is immediate from Christoffel-Darboux formula (Theorem 12.6( 11). 
If Q K {a) = for all k, \k\ = k, we have 



(a (a) - a{y))v = - ^ d 2K A k [n, fi]P K (a)P^(y) G U k+1 
and, therefrom, 



\n\=k 
\fi\=k+l 



a{a)v = T k (v). 

We now show how to obtain a bound for the size of 5-codes, as a function 
of 5. Therefore, in order to cope with any possible 6, we must perturb the 
endomorphism T k as explained next: 
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Theorem 4.4. Let e G M. nk , with e K > 0. Let Tj: be the endomorphism 
defined on Sk by 

T£(v)=T k (v)-e*v k 
where e*v k := J2\ K \=k ^v k Pk- 

(i) T k has a unique maximal eigenvalue At, of multiplicity one, pos- 
sessing an eigenvector v e with positive coefficients. Moreover, if 

Afc-i < X € k < Xk 

(ii) Let e / 0. Any S-code C such that s = m — 5 2 < X e k satisfies 



where a K := Q K (l, ...,!) = T,\n\=k+i A k[ K ,fA- 



Proof, (i) The matrix J k of Tl is equal to except the diagonal elements 
lying in B k . Replacing J k by J k + Mid for some appropriate M, we obtain 
a non-negative matrix which is irreducible so its largest eigenvalue has mul- 
tiplicity one and has an associated eigenvector with positive coordinates. It 
remains true for J k . Since, when e / 0, Jf, < we have A| < Xk- The 
proof of the inequality Xk-i < A| is the same as the one of A^-i < A&. 

(ii) We have av € = X € k v e + e*v\ + u where u G Ilfc+i. We need to compute 
u, and we set u = Yl\^,\=k+i Let ^> I A* I = & + 1> we have: 

\n\<k 
\ K \=k 



and we obtain = Yl\ K \=k v K A k[^-, (A- We have found 
u = Yl ( Yl v * A k[Kil4) p n = Y 

\fj,\=k+l \n\=k |re|=fc 

hence the "generalized Christoffel-Darboux formula": 



(8) ^ e = E w 

s=0 U 



^ - At 
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Now we proceed like in Proposition 12.81 Let the numerator of the right 
hand side be denoted by N k +i(y), and let 

F2 k+ i(y) :=-^-^=JVi H . 1 (y)i;. 

We have: 

f = [F 2k+1 ,l] = [N k+1 ,v t ] 

i«i=fc [/ti=fe 

|re|=fe 

Since the coefficients of e and of v e are non-negative numbers, and /o / 
when e 0, it follows that i^fc+i satisfies the condition (i) of Proposition 
12,31 Condition (ii) is clearly fulfilled if s < Ai. We calculate 

„ M ^ (E| K |=fc<(e K + a«)) 2 

m - A| 

hence the announced bound. 

Let us show that we have indeed generalized the situation described in 
subsection 12.31 and Proposition 12.81 Let x £ M m such that P K {x) > 0, 
Q«(a?) < for all \k\ = k, and P K (x) > for all |«| < k. Let e G R' rfe be 
defined by: e K = —Q K (x)/P K (x). We can show that X e k = a(x). Indeed, 
from (i) and (ii) of the proposition, 

„ . Ei„M<(-W«W/^«W)fl.(») + Q«(»)) 

" {y) = • 

When we let y tend to x, the numerator tends to 0. Since the coordinates 
of v e are positive and P K (x) > for all |k| < k, the left hand side cannot be 
equal to zero when y = x (Pq = 1). So the denominator also tends to zero, 
and A^, = a(x). The Christoffel-Darboux formula (Theorem I2.6f i)) shows 
that v e = T.\ K \<k dl 2 K P K (x)P K . 

When m = 1, n k = 1 and any e > is of this form. When m > 2, it is not 
clear.. It is not even clear that at least one x satisfying these inequalities 
exists. 

Another natural question concerns the values that A| takes. It is hoped 
of course that all values in the interval ]Afc_i,Afc] are attained. We have 
defined a mapping from [0, +oo[ 7r ' s to ]Afc_i,Afc], sending e to A|, which is 
continuous, hence the image in an interval, containing X k , since clearly it 
is the image of e = 0. Let us prove that X £ k tends to A^_i when e tends to 
+oo. To that end, we use the following inequality, valid for any non-negative 
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matrix J with maximal eigenvalue A (|22j): 

(xJ)i 

For all x,Xi > 0, A < sup . 

i 3Ci 

This inequality remains true for the matrix J|, although it is not non- 
negative, because we can apply it to some Jf, + Mid, an argument that 
we have already called for. We choose for x £ M Sk a vector, which first s^-i 
coefficients constitute a positive eigenvector of Jk-i for the eigenvalue \k-i- 
Its last 7Tfc coordinates are denoted by u = {u K )\ K \ =k . We have: 



If \v\ < k 
If \u\ = k 
If |k| = k, 



-(B}.[k,k\ 



The last equality relies on a result that is only proved in Section 5, Propo- 
sition namely that B^[ ] = when k ^ k 1 . 

Let us now choose an arbitrary small a > 0; we can choose the coefficients 
Y.\ K \=k k_ > 1 a ^ v w ith \v\ = k — 1. Then we 



u K > such that 

can choose e K > such that 
are left with: 



Y,\u\=k-i x » A k-iW, K ] 



+ (B k [n, k] - e h 



0. We 



f If \v\ <k-2,^- 
If \u\ = k 



\fc-l 



1,^< Afc_l + Q 



If \k\ = k, 



0. 



Hence A| < A^_i + a for that choice of e. 

Let us go back to the bound proved in Theorem 14,41 We can simplify 
further this bound, getting rid of the eigenvector. We obtain the following 
nicer, but weaker version: 

Corollary 4.5. Let C be a 5-code such that m — 5 2 < Xk—i- Then, 



0) 



\C\< 



4 E| K |=fc d 2KaK 
m - A fc 



Proof. If C satisfies S 2 > m — Afe_i, since Afc_i < A| for all non-negative e 
(from Theorem 14.41 (i)), the bound of Theorem 14.41 (iv) applies to C. We 
get, using Cauchy-Schwartz inequality, and A| < A^: 
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, c] < 1 (S|«|=fc + 



'M 

: ^ d<2K 

\n\=k 



^ _ (e K + a K ) 2 

m-X k 



The function z — > ( z+ct ) is minimized over ]0, +oo[ when z = a. We 
obtain, with e K = a K , the announced bound. 



5. Asymptotic behavior of the largest eigenvalue of T k 

In this section, we compute the limit taken by when the quotient n/k 
tends to some fixed value (Theorem 15.3(1 . This result is needed to pass to 
the asymptotic in the inequality © for the size of a Grassmannian code. 

We first need some very explicit formulas for the coefficients of the sym- 
metric matrix J'& associated to the endomorphism Tk, in the orthonormal 
basis 

{\/d2 K P K , \n\ < k}. From now on we change our usual convention: if not 
specified, tz is a partition of degree s. The diagonal coefficients of J\ are 
the same as the ones of J}., while the other coefficients, denoted by A' s [k, fj], 
satisfy 



A' s [k,h] = A s [n,fi]< 



2 k 



2/i 



To start with, we gather some known results on the polynomials C K . 



5.1. Review of some properties of the polynomials C K . The coeffi- 
cients [£] and 0~) are defined respectively by the following properties: 



oC K = 



H=s+i 



(10) 
(11) 

\u\=s-l 

and have the following explicit expressions: 



(12) 
(13) 



r.(0- 



n 



2ni — 2nj + j — i + 1 
2K/^ 2hij — (— ? % 

i , -, , m ~ i\ tj — 2k j + j — i + 1 
(Ki + 1 H ) 
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while any other values are equal to zero (see |171 Lemma 7.5.7], |14| Th 
14.1], [TH] ). The polynomials C K are intimately related to the decomposition 
of GL(m,R)-modules (0 Theorem 5.2.9]): 

ft(GL(m,R)/0(m,R)) = 

For later use, we settle the notation: 5 K := dim(i 7 ^) and we recall the 
formula (0): 

Ki Kj + J I 



(14) 



S K := dim(F^; 



n 

l<i<7<m 



5.2. Formulas for the coefficients of the matrix J'k- 

Proposition 5.1. The matrix B s has the following properties: 

(i) B s [k, k'] = for all k ^ k' . 

(ii) Ifm < n/2, 



2B s [k, k] = 



o 




2k{ + m + 1 — i 




K 


2Ki + n/2 + 1 - i 



(hi) If m = n/2, B s [k,k] = m/2. 





K 


2ni + m — 1 — i 


\ K (i)J 




2k { + n/2 -1-i 



Proof. We recall that the coefficients fi K ^ u are defined by: 

v\k>v 

Inverting these relations, we obtain coefficients a KV such that 

C R = a K P K + ^2 °k,uPu- 

v\k>v 

Taking into account the formula (|1U|) . we obtain: 



B s [k,k'} = K ^2 ^ a W + E @ K 



Oi K '- 



\fi\=s+l \u\=s— 1 

We use the following obvious relations: ol k (5 k — 1 and oi K f5 K ^ K i + Pk>ol kk i 
to rewrite 



(15) 
(16) 



B s [k,k'] 



E 

IH=*+i 



Pu,K,' 

Pa 



+ E 

\u\=s-l 



Pn, 



K Pn 
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Let us assume first that k ^ k'. Since rz\ is non zero only if ji 



for some index i, and also 0uk' is non zero only if fj, = k'^' for some index 
j, at most one term in the first summation may be non zero, and the same 
argument holds for the second summation. We only have to consider the case 
when k' satisfies: for some indexes i 7^ j, k[ = Ki + 1 and = Kj — 1. The 



0') 



remaining terms in the expression of B s [n, k'} correspond to fi 
and v = K(j\ = k' 'm. 

Moreover, the coefficients [5^ )K are calculated in |14| . and in particular we 
have: 



(17) 



K 

K U) 



2k j + m — 1 — j 
2nj + n/2 -I- j' 



Replacing in (|15|) we have 

A, 



B s [k, re'] 



2Kj + m — 1 — j 



P K > 2(2 Kj + n/2 - 1 - j) ■ 



.(<) 



Combining (|1U|) and Qll[) in the obvious relation: 

(ecr - cre)C K = mC K 



leads to: 
(18) 



E 



E 

|j/|=a-l 













-i 















if re / k' 
m if re = re' 



From l|18|) we can conclude that B s [k, re'] = when re 7^ re'. 
When k = k', replacing ((17)) in l|15|) leads to the formula (ii). If moreover 
n = m/2, taking account of (|18|) we obtain 5 s [re, re] = m/2. 
We now give explicit formulas for the coefficients of «/'&: 

Proposition 5.2. With the following notations: 

qi := 2Ki — i + m 



D(x) 



N := n — 2m 

x 2 



x 2 -l 

(x+l)(x+N) 



and 



U W ~~ (2a:+ AT) [2x+N+2) X T 1 u 

(7(0)- 1 



7V+2 
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we have the expressions: 
B s [k,k] = 



t-t £(LT^.-* + i)); * + 2 



iV 



2 qi + N + 2 



ied(«) i=i 



2ft + iV - 2 



A'. [«,««] = ((]jD(ft-ft + l)L>(ft+ft+iV + l)). 



Proof. For the calculation of B s , we replace in Proposition 15.11 (ii) the for- 
mulas (j!2|) and (jL-jj) , and take account of Proposition 15.11 (iii) . 
In order to calculate A' s [k, k®], we have already seen that: 



A s [k,k® 
We need a formula for 



-a — i • Expressions for the leading coefficients 

of the polynomials C K and P K can be found in ^2] and |2Sj. Putting them 
together we find: 

_ JT / ft + ft + iV \ (2ft + AQ(2ft + iV + 2) 

where the last fraction must be understood as (N + 2)/ (N + 1) when ft = 0. 
Joined with (1101). we obtain 



j4 s [k, k 



01 



n 



% - ft + 1 \ / ft + ft + jv + 1 



j=i \ ii- ft 



Next we use (0): 



2k 



n 

3^1 



ft + ft + N 

(q l+ N)( qi + N+l) 
(2ft+iV)(2ft + 7V + 2)' 

ft-ft + 2\ /ft + ft + 7V + 2 



Qi ~ Qj J V % + ft + AT 

(2ft + jV + 4)(ft + AQ(ft + iV + l) 
(2ft + A0(ft + l)(ft + 2) 
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where the last fraction must be understood as (A r +4)(A r + 2)/2 when qi = 0, 
and we obtain the announced formula for A'[k, k®], 

(n) 

5.3. The limit of A/%. Now n varies with k so we rather denote by the 
endomorphism defined previously and A^ its largest eigenvalue. 

Theorem 5.3. If n/2k — » I, while n — > +oo and k — * +oo, 

A („) _ J+l/m_ 

Proof. We give careful proofs in the cases m = 1 and m = 2, and will be 
more sketchy in the general case. As it was noticed previously when m = 1 
the eigenvalues are the zeros of the Jacobi polynomials; their asymptotic 
is calculated in [7j, exploiting the differential equation for the Jacobi poly- 
nomials and Sturm's method. Another approach, using chain sequences, is 
used in [21]. However, none of these methods seem to generalize easily to the 
several variable case. Our argument will only use the fact that the matrix 
J'^ is non-negative. More precisely, we use the following: 

Lemma 5.4. [22] Let J be a non-negative symmetric matrix of size N, with 
largest eigenvalue A. 

(i) For all x G R N with xi > 0, A < max; 

(ii) For all x £ R N x + 0, A > <£ ^ L . 

The case m = 1. We recover from Proposition 15.21 the formulas: 

(n - 2)(n - 4) 



26, = 1 



(4s + n)(4s + 7i — 4) 
(2s + l)(2s + 2)(2a + n - 2) (2a + n - 1) ^ 1/2 



(4s + n - 2) (4s + n) 2 (4s + n + 2) 
From these expressions we see that both sequences are increasing with 
s. Moreover, we see easily that if s ~ k, b s ~ 2 rapgp , and a' s ~ JJ+Tj?- 
Applying Lemma 15.41 (i) with x s = 1 for all s leads to: 

Ajf ) < a' fc -i + 6fe + a'k 
and the right hand side tends to 4(£ + l)/(£ + 2) 2 when n/2k tends to i. 

We lower bound A^ using Lemma 15.41 (ii) and a choice of x proposed in 
: let x be defined by: 

x s =0 1 < s <t:=k- [Vk\ + 1 
x s =1 t + l<s<A; + l 

so that x s = 1 on the [v^J last coordinates. Then, 
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(xJV)-x > T,:Z} + l(*'s-l + b s + a> s ) 



x • x k — t + 1 



Again, the right hand side tends to 4(£ + l)/(£ + 2) 2 , hence the result. 
The case m = 2. From Proposition 15.21 we have, setting s := n\ + K2 and 

V := Ki — K2'. 

(n-6)(n-4) { (4s + 2n-4) 



8(4s + 2n - 6) V (4«i + n)(4K 2 + n - 2) 
(4s + 2n - 8) 2 
(4ki + n- 4)(4k 2 + n-6) 

A s [k,k 



(Dl _ f_(2^+2)M 1/2 / (2.s + n-2) 2 ^ ' 2 



(2v + 2) 2 -lJ \(2s + n-2) 2 

(2/ci + 2)(2ki + 3)(2ki + n - 3)(2k x + n - 2)y /2 
(4ki + n - 2)(4ki + 7j) 2 (4ki + ra + 2) 

As[K ' K l ~{(2v) 2 -l) {(2s + n-2) 2 - 

(2k 2 + 1) (2« 2 + 2) (2« 2 + n- 4) (2k 2 + n - 3) \ 1/2 
(4K 2 + n-4)(4K 2 + n-2) 2 (4K 2 + n) J 

One can verify that these coefficients are increasing with s when v stays 
constant. This is easy to see for B s , not so obvious for the two others because 
the second term is decreasing while the last big quotient is increasing. 

In order to obtain a lower bound for from Lemma 15.41 (ii). we choose 
x = (x K ) with: x K = 0, 1. We fix a number V < k — \fk. Let /Cy,s be the set 
of the V partitions of degree s with smallest v = K\ — k%. Hence ICv,s = { K I 
\k\ = s,k 2 > Lf J - V + 1}. We set x K = 1 iff deg(re) > t := k - [v^J + 1, 
and k £ /Cyui. We need to avoid in K\r,s some partitions, namely the ones 
with v = and the ones with v maximal (for those partitions, some terms 
are either missing or are equal to zero in (xJ'^) K ). Let this new set be 
denoted by JC'v s - We have, when x K = 1, k £E K'vsi N k, 



i=l 
2 

i=l 
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In the expressions of A' s [k, k®] we can minor the first term by 1 (v 7^ 0), 
then minor each term by its minimal value in the sequence v = cte to which 
it belongs. As was mentioned before, this minimal value is obtained when 
the degree is minimal, i.e. when s = t or s = t + l. We do the same for 
A' s _x[ K (i)) K ] an d for B s [k, k]. Then we must consider the behavior when s 
is constant of B s [n, n], of: 



((2a + n- 2) 2 - 1)1/2 
(C( qi )C( gi + l)) 1 / 2 + (C(g 2 )C(g 2 + I)) 1 ' 2 



and of the analogous expression .A 2 s _i[k] corresponding to the last term. 
These expressions are increasing with k 2 . Let B s y mtn , A^ s y be their 
minimal values in IC'vs- F° r simplicity, we assume that 

™\-^(T3 mm R min\ rj min 

mm{B t ,V , Bt+i,V ) — B t y , 

and the same for A 1 , A 2 . 
We obtain: 



x ■ x 

(k - t - l)(V — 2) min, 4 1 min . 2 mirl \ 

(k-t + i)v {Bty +Aty +Aty ] - 

Now we let n/2k tend to I. Since Bty mm is obtained at a partition 
essentially equal to [t/2 — V/2,t/2 + V/2, and since i ~ k, we see that 
£ f ,y min tends to 2(1 + l/2)/(i + l) 2 . For the same reason, A x t y min and 
^mm tend tQ ^ + 1 / 2 )/(£ + 1)2 (the parameter V is still fixed at this 
stage). So we obtain 

liminfA^>(l-|).4^. 

Now we let V tend to +00 to obtain the appropriate lower bound. 

The second and last step obtains an upper bound for Aj^ from Lemma 
15.41 (i) with an appropriate choice of x. The choice x K = 1 for all k is not 
good enough here because D(2v + 2) 1 / 2 + D{2v) x l 2 7^ 2. We need some 
x K that modify properly these factors. We choose x K := (2v + l) 1 / 2 where 
v = K\ — K2- We have 
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( T/(™)\ 2 

= B a [ K , k] + V K »]5s« 



i=l 



Let: 



( 2v + 2 \ f (2s + n-2) 2 \ 1/2 



1/2 



Since 



/(2ki + 2)(2ki + 3)(2«i + n- 3)(2ki + n - 2) 
\ (4«i + n - 2) (4/ci + n) 2 (4«i + n + 2) 

' ( 2« 2 + l)(2/c 2 + 2)(2k 2 + n - 4)(2« 2 + n - 3) \ 1/2 
(4k 2 + n - 4) (4k 2 + n - 2)2 (4k 2 + n) ) 

A>\ 1 f 2v U (2a + n-4) 2 \ 1/2 
= l^TlJ ( v (2,s + n-4) 2 -lJ • 

/(2ki)(2ki + l)(2/ti + n - 5)(2«i + n - 4)\ 1/2 
\(4ki + n - 6)(4«i + n - 4) 2 (4ki + n-2)J 

. br . /2« + 2\ / (2 S + n-4) 2 \ 1/2 

( (2« 2 - l)(2/c 2 )(2/c 2 + n - 6)(2/c 2 + n - 5) \ 1/2 
V (4k 2 + n - 8)(4k 2 + n - 6) 2 (4k 2 + n - 4) / 



^ + ^ ( 2v + l)V2 " 2v + 1 > 



and 



we have: 



2v 



(2w + l) 1 /2 2v + l' 



fc 



/ j/\ n )\ 



i=i 

2 

+ J]A b [K (i) ,«]. 

j=l 

This expression is increasing with s when v is fixed. When s is fixed, 
Bs[k,k], £- = iA»[/e,K«] and E?=i Ab [«« 

, n] are maximal at k = [s/2, s/2] 
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(we extend the functions to partitions with real parts here). We obtain, 
with Pk = [fc/2,fc/2], 



_k 

K x* 



max - < 



B s [pk,Pk] +Y; Ai lPk,Pk {i) ] + ^2^[p Hi) ,p k ]. 



i=i i=i 



The computation of these values shows that the right hand side tends to 
A(l + 1/2) /{i + l) 2 when n/2k -► 1 

The general case m > 2 works the same. For the lower bound, we use 
K-V,s '■= {k I \k\ = s,K m > |_^J — V + 1}. The cardinality of Ky, s only 
depends on s mod m. We should avoid some partitions in /Cy iS , namely the 
ones with some parts equal and the ones with n m = [^J —V + l. Their 
number is negligible compared to the cardinality of /Cv>. Then, we proceed 
in the same way as for m = 2. 

The upper bound is obtained with x K = (52k) 1 / 2 . We have 

m / c \ 1/2 m . -, 



7=i v ° 2k / j=i « " n 



hence 



Xr i=i 

m 

+ ^2 A h [k^, k], 



i=i 



where 



A*[«,«W] = fT*Z^±l. 
fJi U " '// 



m 

[] + + AT + l))C(q t )C(q t + 1)) 



1/2 



J'=l 
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and similarly 



11 Qi-Qj 



II D ^ + 9j + N~ l))C( qi - 2)C( qi - 1) ' ' 



We have the nice identity: 

- gj i + 1 



En 



m. 



, % - Qj 



We do not have a reference for this last identity, so we give an argument 
here: from (|14j). 

dimF^ _ -j-r - Kj + j - i + 1 
dim i 7 ^ -fj- Kj — Kj + j — i 

We obtain the demanded identity as the equality of the dimensions in the 
following decomposition of GL(m, M)-modules (Pieri's rule, [S]): 

f$®fz 1 = ®™ 1 f£ ) . 

It turns out that the coefficients B[k,k], A$[k, «W] and A 1, [rem , k] are 
increasing when k runs over a sequence of the type {y + 1, . . . , l]) s >o 
(when N is big enough), and that, on the space of partitions (with real parts) 
k of fixed degree k, the maximum of the expressions B[k, k], Y^lLi A$[k, kW] 
and YliLi j 4 t, [ K (i)i K ] is attained at k = = [k/m, k/m, . . . , k/m]. 

Moreover, it is easy to see that, when n/2k — » £, 

v r r 1 o *+V"» 



+ 2/m) 2 

(m \ / m N 

2><WS?] = lim J^W^ft] 
i=l J \i=l > 

£ + 1/m 
~ (£ + 2/m) 2 ' 



Remark 5.5. One obvious consequence of Theorem \5.cl\ is that, for fixed n, 
the eigenvalue Al n ^ runs over the whole interval ]0, m[. Hence the bounds 
proved in section^ for the size of grassmannian codes potentially cover all 
possible minimal distance. 
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6. An asymptotic bound for the size of Grassmannian codes 

We are now ready to take the limit when n tends to +00 in the inequality 
Q, and prove Theorem ll.il 

We are left with the estimate of logQ^ui_ fc d2 K a n)/ n - 

Lemma 6.1. Let 5 R := dim(i^). Ifn/2k — > S M while n and k tend to 
+00, 



(19) limsup— log( ^ 5 K ) < 



n 

\n\=2k 
e(n)<rn 

(20) m((l + p)log(l + p)-plog(p)). 

Proof. In the case m = 1, 62k = dimS^fc = ( n+ 2fc _1 ) anc ^ ^ ^ s a classical 
result. The general case is probably well-known but since we lack a reference, 
we give a proof here. Let k be a partition of length at most m and of degree 
2k, that we extend to a partition with n parts with an appropriate number 
of zeros. From 1)14)1. 

dim(^)= n 

l<i 

Since Kj = when j > m, we have 



7 — z 

l<i<j<n J 



Kj Kj + J 2 -1— r "I T + J — I 



n ^"7-7" n n 

We upper bound: 



l<i<j<m l<i<m j>m 



n 



«i + j - i fn + Ki-1 



J — 1 



and 



to obtain 



TT Ki Kj+j i <(2k + i) 

A ± 7 — i 

l<«<i<m 



|re|=2fc |re|=2fc \i=l V * 7 

e(K)<m £{K)<m 

2k 



^( 2fc+i r 2 (E( n+ r 1 ) 

< (2fc + 1)' 



s=0 

n + 2k 
2k 
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and we obtain the announced limiting result using the classical 



1 [n + 2k\ . , 
hm -log = (l + p log (1 + p -plog(p . 

n/2k^p- 1 n \ Ik J 

From the three-term relation (|2.4|) . specializing to (1, . . . , 1) we get triv- 
ially 

a K <m and hence Yl\ K \=k d2 K a K <mE H=2fc (obviously d 2/t < <5 2k since 

£(n)<m 

V£ K is contained in F^ K ). 

Then we only have to solve the equation, involving the limiting result of 
Theorem 15. Ml 



lim A 



(n) 
fc-1 



which leads to 



m , 

P=y(-1 + (1 



p 1 + l/m 

- -r 1/2 )- 

m 



7. LP versus Hamming 

In 4 , A. Barg and D. No gin give an asymptotic bound for the size of 
Grassmannian codes, derived from the so-called Hamming bound. They 
prove, with the notations of Theorem ll.il 



Theorem 7.1. [3j 

(21) 



ilog|C|<- m .o E ( 1 /l- v /l±p 



It turns out that our bound is better than (|21|) only when s is small. 1 
The crossing point sq for the two bounds has the approximate value: 



m 


2 


M 


4 


5 


6 


so 


1.4528 


1.2714 


1.1853 


1.1372 


1.1067 



7 


8 


9 


10 


1.0856 


1.0702 


1.0584 


1.0492 



Figure n plots the two bounds for m = 2. 
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Figure 1. LP and Hamming asymptotic bounds for m = 2 
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